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ABSTRACT 

We constrain the ratio of black hole (BH) mass to total stellar mass of type-1 AGN in the COSMOS 
survey at 1 < z < 2. For 10 AGN at mean redshift z ~ 1.4 with both HST/ACS and HST/NICMOS 
imaging data we are able to compute the total stellar mass M^totai, based on restframe UV-to-optical 
host galaxy colors which constrain mass-to-light ratios. All objects have virial MBH-estimates available 
from the COSMOS Magellan/IMACS and zCOSMOS surveys. We find within errors zero difference 



between the M-qjj-M.. 



total' 



relation at 



1.4 and the M BH -M : 



, bulge' 



relation in the local Universe. 



Our interpretation is: (a) If our objects were purely bulge-dominated, the Mbh _ M«, bulge- relation has 
not evolved since z ~ 1.4. However, (b) since we have evidence for substantial disk components, the 
bulges of massive galaxies (M*. total = 11.1 ± 0.3 or logMeH ~ 8.3 ± 0.2) must have grown over the 
last 9 Gyrs predominantly by redistribution of disk- into bulge-mass. Since all necessary stellar mass 
exists in the galaxy at z = 1.4, no star- formation or addition of external stellar material is required, 
only a redistribution e.g. induced by minor and major merging or through disk instabilities. Merging, 
in addition to redistributing mass in the galaxy, will add both BH and stellar/bulge mass, but does 
not change the overall final Mbk/M*, bulge ratio. 

Since the overall cosmic stellar and BH mass buildup trace each other tightly over time, our scenario 
of bulge-formation in massive galaxies is independent of any strong BH-feedback and means that the 
mechanism coupling BH and bulge mass until the present is very indirect. 

Subject headings: galaxies: active — galaxies: nuclei — galaxies: evolution — galaxies: fundamental 
parameters 



1. INTRODUCTION 
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Masses of galactic bulges and their central black holes 
(BHs) follow a tight re lation in the local Universe (e.g. 
iMarconi fc Huntl 120031 : lHaring & Rixl 12004 ) with only 
0.3 dex scatter - strong evidence for a coupled forma- 
tion and evolution of galaxies and BHs. The source of 
this coupling is unclear, but feedback mechanisms have 
been proposed involving the central potential well depth 
regulating BH accretion, or more violent feedback from 
active galactic nucle i (AGN) into their host galaxies (e 
[ Hopkins et~al]|2006t ISomerville et ail 12008 PMenci et al 
l2008f) . While these scenarios potentially provide ingredi- 
ents for acquiring consensus with observations, all such 
models include ad hoc assumptions and do not work 
from first principles. Empirical constraints are urgently 
needed to investigate the actual physical processes in- 
volved in the coupled evolution. 

One strong constraint is the evolution of the Mbh~ 
-W*, bulge-relation over time. While circumstantial evi- 
dence grows that the va lue of Mb^/ M^^ ^p. was larger 
at earlier cosmic times (|Peng et al1f 2006a.b; Treu et al 
[21)071 [Woo et al.|[2008l:IWalter et al.ll2004t IRiechers et al 
2008a, b, 2009), studies are subject to biases (jLauer et al 
120071 ) and better statistics are required to investigate 
where in Mbh, or when in cosmic time, a turnoff from 
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the local Mbh - -^*, bulge-relation occurs. 

Broad-line AGN and their host galaxies are the only 
systems at higher redshifts in which both the mass of 
the galaxy or its bulge as well as its central BH can 
be estimated. Here we set constraints on an evolving 
Mbh-M* bulge-relation by computing optical color-based 
stellar masses (from HST/ACS and HST/NICMOS) and 
combine them with virial Mbh (from Magellan/IMACS 
and zCOSMOS/VLT/VIMOS) for 10 AGN in the red- 
shift interval 1.06 < z < 1.92, 3.2-5.5 Gyrs after the Big 
Bang. 

Throughout we use AB zero-points unless otherwise 
noted and a cosmology of Hq = 70kms _1 Mpc _1 , = 
0.3, and Q A = 0.7. 

2. DATABASE 

In order to control selection effects on M*, Mbh or 
any special relation between these two, we require a 
transparent sample definition. Our selection of type-1 
AGN is based on X-ray detection in the XMM-COS MQS 
survey (jHasinger et al.ll2007HCappelluti et al.ll2009T ) and 
subsequent identifi cation of their optical counterparts 
(jBrusa et al.l 120071) . Classification as type-1 AGN for 
this study uses both spectroscopic identific ation of broad 
emissi on lines in the M agellan/IMACS (iTrump et al.l 
[2009ah and zCOSMOS (| Lilly et alJ 12003) surveys as 
well as photometric classification usi ng the long spectral 
baseline SEP cove r ed in COSMOS (ICapak et al.l [20071: 
ISalvato et all 120091 ; Elbert et al.H2009al Capak et al.. in 
prep.). As such, only objects with high-confidence clas- 
sification enter our sample. 

For ~550 type-1 AGN selected this way we require 
the following data, resulting in a random subsam- 
ple: (1) coverage by A CS F814W (IScoville et al. 1 120071 ; 
IKoekemoer eTall[2007l) . and (2) NICMOS parallels 16 , as 
well as (3) spectra for virial MeH-estimates from the 
Mgll broad emission line. 

The limiting factors are the relative NIC3 coverage of 
6.4% of the ACS area (~20 AGN) and the status of ongo- 
ing spectrosc opic surveys. Black h ole masses have been 
calcula ted bvlTrump et al.l (2009bD from IMACS spectra 
and by IMerloni et al.l (|2009f) using zCOSMOS/VIMOS. 
For this letter we have BH masses available for 10 AGN, 
spanning the redshift range 1 < z < 2. Two AGN were 
observed by both instruments - the MeH-estimates are 
consistent within 0.2 dex in both cases. Sample informa- 
tion is listed in Table [1] including BH masses and galaxy 
parameters derived below. 

3. HOST GALAXY MASSES 
3.1. Observed host galaxy photometry and colors 

We obtain information on the host galaxies using 
broad-band photometry from the high-resolution HST 
ACS/WFC images i n the F814W (=1) fi lter with 
O'.'03/pixel sampling (Koekcmo er et ail [2007). and the 
NICMOS/NIC3 parallels in the F160W (=H) filter at 
O'.'101/pixel, both integrated for one orbit. 

In the NICMOS i?-band the host galaxies of all AGN 
are clearly resolved, visible already to the unaided eye. 
We hence extract the host galaxy flux from the compos- 
ite galaxy+AGN NICMOS image by modelling the two- 

16 http:/ /irsa. ipac.caltech.edu/data/COSMOS/imagcs/nicmos/ 



dime nsional light distr i bution of e ach object using GAL- 
FIT (|Peng et al.H200l IPendl2009T) in Version 3.0 (C. Y. 
Peng, priv. comm.). We restrict ou r mod els to a point- 
source plus a single elliptical Sersic (1968) profile. Previ- 
ous simulations show that at our resolution and depth it 
is unre liable to use many and/ or complex galaxy c ompo- 
nents (Simmons fc Urrvl [20081; ISanchez et al.ll200l . We 
carry out several passes of GALFIT, first with free Sersic 
parameter 17 n, and subsequently fixed n = 1, 2 or 4, de- 
pending on the best free-n-fit. 

We require GALFIT to converge on a sensible solution, 
indicating that an actual host galaxy is being described 
and not e.g. uncertainties in the point spread function 
(PSF). This means we require a GALFIT nucleus-to-host 
contrast of < 3.25 mag, a half-light radius of r > 2 pixel, 
and 0.5 < n < 8. 

XID 2261 and 5230 both show unrealistically compact 
host galaxy model scale-lengths, indicating an unsuccess- 
ful GALFIT model. For these two cases we simply sub- 
tract the best fitting single point-source model - without 
a Sersic component - from the original image, resulting 
in an only slightly oversubtracted/underestimated host 
galaxy. We use aperture photometry on the host galaxy 
in these cases, and the GALFIT galaxy model magnitude 
for the other eight. 

In total we resolve all 10 host galaxies in the sample 
- no object drops from the sample due to high nucleus- 
to-host contrast or other reasons. Extracted host galaxy 
images are shown in Figures |T]+[2] 

In the ACS J-band data the contrast between AGN 
and host galaxy is less favorable than in the i?-band, 
as expected from the near-UV SED. Therefore we take 
a two-step approach. First, we carry out "peak sub- 
traction" removal of the nuclear component for the I- 
band, by scaling a PSF to the central 4 pixel aperture 
flux. From statistics on the expected random residuals 
after subtraction of t his scaled PSF for several 1000 stars 
(| Jahnke et al.ll2004bl ). we require a residual flux of >5% 
for a host galaxy to have a high probability of being re- 
solved. This is the case for 7/10 objects - the hosts of 
XIDs 59, 2261, and 5230 remain unresolved in the ACS 
image. 

The seven resolved objects are again modeled using 
GALFIT, with Sersic n fixed from the iJ-band fit, thus 
minimizing S/N-dependent biases in our extracted col- 
ors. As for the iJ-band we check the models for success- 
ful convergence, which is the case for all seven resolved 
objects, and again use the host model magnitude. The 
resulting host galaxy photometry is also listed in TableHJ 

3.2. Host galaxy stellar masses 

Our general approach is to constrain the mass-to-light 
ratio (M/L) of each host galaxy by a single optical color. 
We have successfully employed this method before to 
study stellar popula tions of low-z quasar host galaxies 
(jJahnke et al.ll2004al ) or stellar ages and masses of quasar 
hosts a t 0.5 < z < 3 (ISanchez et al.ll2004t Uahnke et al.l 
I2004ht IKuhlbrodt et all 120051 : ISchramm et alJ 12008ft . 
Here we compute stellar masses from the rest-frame V- 
band luminosity in combination with the M/L from the 
rest-frame (B — V) color: 

17 n = 1 represents an exponential disk, n = 4 a de Vaucouleurs 
spheroid. 
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TABLE 1 

Sample summary 



XMM-Newton name" 


XID a 


z 


M BH b 

(Mq) 


Ref 


F814W host 
(AB mag) 


F160Whost 
(AB mag) 


Sersic ra d 
F160W 


(B - vy%* 

(Vega mag) 


log(M„ itotal ) e 
(M ) 


XMMU J100118. 5+022739 


14 


1.065 


8.52 


1 


20.70 


19.41 


1.9 


0.38 


11.18 


XMMU J100046. 8+020016 


59 


1.923 


8.72 


2 




20.28 


0.6 


0.15-0.75 


11.07-12.12 


XMMU J095927.7+020010 


219 


1.248 


8.07 


2 


22.55 


20.61 


2.1 


0.49 


11.00 


XMMU J100035. 3+024303 


281 


1.177 


8.25 


1.2 


22.63 


20.25 


3.5 


0.74 


11.44 


XMMU J095928. 5+015934 


329 


1.166 


8.05 


1,2 


22.95 


21.03 


1.7 


0.58 


10.90 


XMMU J100130.7+021147 


2148 


1.526 


8.43 


1 


23.45 


21.15 


1.5 


0.40 


10.88 


XMMU J100243. 8+020502 


2261 


1.260 


8.05 


1 




21.38 




0.15-0.75 


10.25-11.09 


XMMU J100226. 9+015938 


2637 


1.630 


8.35 


1 


24.12 


20.72 


1.8 


0.73 


11.64 


XMMU J095903. 2+022001 


5049 


1.131 


8.40 


2 


22.66 


20.82 


1.5 


0.57 


10.93 


XMMU J095908. 1+024310 


5230 


1.359 


8.22 


1 




19.13 




0.15-0.75 


11.21-12.08 



Original XMM-Newton source name and ID (Cappcllut i et al.|[20Qgh 



Mean value where two measurements are available 



bourcc for M BH : (1) Magcllan/IMACS (Trump ct al. 2009b); 
0.4 dex and 0.3 dex, respectively 

From free-n fit before fixing 
c Total uncertainty is ±0.4 dex (Scction[3?3j . 



(2) zCOSMOS/VIMOS JMerloni et~al| [2009'); M BH -crrors arc quoted as 



XID 14 



XID 59 




XID 219 


XID 281 

• 

• 


XID 2637 

• 


XID 5049 

1 



Fig. 1. — Nucleus-removed host galaxies: Galaxy plus nucleus model fitted for 8/10 objects (HST/NIC3 F160W). Images are 7" X 7", 
some objects lie near NICMOS tile edges. 



XID 2261 fc XID 5230 

* ft 



Fig. 2.— As Figured] XID 2261 and 5230, only nucleus model 
fitted and subtracted. 

L V = io-°- 4 ^- 4 - 82 ) (1) 

M*, total = 10-O.M2+1.710(fl-V) x (2) 



with L and M in solar units. This calibration 
is ba sed on template-fitted masses (|Bruzual fc Chariot! 
120031 Chabrier IMF) and lu minosities derived for galax- 
ies in the COSMOS field (|Ilbert et alJ l2009bfl . We 
convert luminosities to Vega zeropoint and apply a 
—0.124 dex mass offset to transform to the mass scale of 
the models by Chariot & Bruzual (2007/2009 in prep.) 
that include contributions from TP-AGB stars. The 
linear relation of eqn. [2] is a fit to galaxies in redshift- 
(1.0 < z < 1.6) and color-range (0.38 < (B - V) < 0.74) 
of our 7 host galaxies resolved in both HST bands. The 
scatter of the fit corresponds to an RMS uncertainty in 
resulting stellar mass of ±0.3 dex. 

We convert our measured (F814W-F160W) colors to 
restframe (B — V^)host by applying both JT-correction 
and a color term, individually for each object and its 
rcdshift. For this purpose we identify for each galaxy 
the single stellar population model again from Chariot 
& Bruzual (Chabrier IMF, solar metallicity) with the 
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closest (F814W-F160W) at a given z. Since the inter- 
polation intervals are rather small and we are not using 
the interpolation SED to extract any further informa- 
tion, this method is quite insensitive to the exact choice 
of models, and errors from M/L-calibration dominate. 

For the three host galaxies only detected in F160W we 
have to make assumptions for the interpolation template 
to convert observed F160W to Ly, as well as the color 
(or M/L) in equation [2] On the red/old side we assume a 
value of (B — V) = 0.75 (Vega ZP) co rresponding to th e 
red end of the red sequence at z ~ 1.5 (jKriek et al.l l2008). 
a value also consistent with the reddest values of the rest 
of the sample. As a blue limit we use (B — V) = 0.15, 
corresponding to 3a from the mean value for COSMOS 
inactive galaxies at z ~ 1.4, logM^totai > 10.7 and age 
t < 1 Gyr. These assumptions should bracket the true 
M/L, and provide robust limits on the host galaxies' stel- 
lar mass. 

3.3. Stellar mass uncertainties 

Uncertainties in stellar mass have several sources. In 
addition to the (B — V) calibration stated above, the 
strongest contributions come from (1) GALFIT precision 
of extracting the host galaxy, (2) potential influence from 
different spatial resolutions in the ACS and NIC3 images, 
and (3) dependency of bandpass conversions on assumed 
temp lates. 

(1) ISanchez et all (pOOl ) tested how reliably GALFIT 
can derive host galaxy photometry in one-orbit ACS 
data. Given the typical host galaxy magnitude of our 
sample in the F814W filter we conclude an uncertainty 
of 0.15 mag for F814W, and - due to counteracting effects 
of lower spatial resolution but more favorable galaxy-to- 
nucleus contrast - also for F160W. 

(2) GALFIT 's functionallity depends on the spatial 
resolution of an image and the spatial difference between 
the AGN and galaxy component. If a galaxy is compact 
and the PSF not well characterized, flux transfer between 
components is possible, its amplitude potentially depend- 
ing on spatial resolution. We test if the different resolu- 
tions of the ACS and NIC3 have a significant influence, 
by rebinning a mock AGN+galaxy image resembling a 
typical object to different spatial resolutions. The recov- 
ered host galaxy photometry shows an rms variation of 
^0.15 mag, no systematic offset, and a negligible trend 
with resolution. This rms scatter is consistent with the 
GALFIT precision from (1) above and we conclude that 
resolution effects are insignificant. 

(3) The bandpass conversion and if -correction de- 
pend on the IMF and metallicity of the assumed sin- 
gle population interpolation model. The masses derived 
from Chabrier and Salpeter IMF generally differ by less 
than 0.1 dex. The metallicity of galaxies of the esti- 
mated masses at 1 < z < 2 is expected to range be- 
tween solar and 3x solar (8.8 < 12 + log( 0/g) < 9.2, 
iTremonti et~all 12004 iMaiolino et~aTl I2008D . Changing 
from the solar metallicity used, to the most metal-rich 
templates (Z = 0.05), masses change by <0.05 dex. We 
conclude that our stellar masses should be good to within 
^0.1 dex from the choice of bandpass conversion SEDs. 

In combination we find our error budget in stellar mass 
is 0.21 mag for the (B — VQ-color (0.15 mag from each 
band), corresponding to 0.27 dex uncertainty in stel- 




Fig. 3.— M B H-M, itotal -relation from COSMOS ACS+NIMOS: 
Shown are ten type-1 AGN in the redshift range 1.06 < z < 1.92, 
seven with direct stellar mass estimates (red circles) and three with 
bracketing range (blue circles and lines). Th e points are overplot - 
ted over the local A/bh~^*, bulge - 

relation by Haring & Rix ( 2004) 
(black triangles) and its best fit (solid line, logMBH = 8.2 + 1.12 X 
(log Af t bulge — 11))- At a given mass there is no difference in the 
A^BH/-^*,total- ra ti° at 2 ~ 1.4 and local MBH/A4*,bulge- ra ti°i f° r 
the sampled BH mass range, logAfaH ~ 8.3. 

lar mass - dominating the three sources of uncertainty 
above. Adding the uncertainty of the mass calibration 
from equation [21 the total uncertainty in stellar mass is 
±0.4 dex. 

4. RESULTS 

While we can not estimate bulge masses directly, we 
find that the seven objects with direct total stellar mass 
estimates fall directly onto the h -M r hi . l fffi - r elat io n 
(Fig. [3]) of the local universe from lHaring fc Paxl (|2004h . 
The three objects with a bracketing range on stellar 
mass are also consistent with the local relation. The 
objects have a maximum deviation of ~0.3 dex per- 
pendicular to the z — relation. The seven non- 
limit AGN in the z-range 1.06 < z < 1.65 show a 
mean ratio M B H/M* jt otai = 0.00178 ± 0.0012 and mean 
log(MBH) = 8.31. This is consistent with the value at 
z = of Mbh/-^*, bulge = 0.00165 at the same Mbh, and 
h as exactly the s ame .3dex scatter. 

iMerloni et al.l ((2009) compute stellar masses for a 
larger sample of COSMOS type-1 AGN, using an in- 
dependent SED-decomposition method. For 18 galax- 
ies where stellar masses could be estimated with both 
methods (five objects are part of this study, 13 have no 
BH mass estimates yet) their masses are smaller by 0.1- 
0.2 dex. This agreement within the error bars reinforces 
our conclusion that our mass estimates are robust. In 
total they find a mild deviation from the z — relation, 
but most of their signal comes from objects at z > 1.5, 
which is not well covered by our study. 

5. DISCUSSION & CONCLUSIONS 
5.1. Completeness and systematics 
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Are there systematic effects inherent in the dataset 
that prevent us from detecting objects which deviate 
more strongly from this relation? X-ray selection and 
subsequent broad-line AGN classification finds all AGN 
down to a contrast of AGN/galaxy<10%, beyond which 
the stellar light swamps any AGN signature. The general 
limit is set by the magnitude limit of the spectral fol- 
lowup, modified by this contrast. At the o bserved Mbh 
we ex pect generally high accretion rates (e.g. lZheng et al.l 
2009) , thus this detection limit converts to a limit in Mbh 
present in the data. The additional ACS and NIC3 cov- 
erage are random processes. Being able to extract a host 
galaxy for every single object of our sample means that 
aside from the MBH-hmit we have no completeness limi- 
tations with respect to galaxy mass. In summary, we do 
not expect missing objects in the MBH-regime currently 
populated. 

We also co mment on the lu minosity function or 
"Lauer" bias ([Lauer et alJ l2007f ): The M B H-Af*, tota i- 
relation has a scatter and the BH luminosity function 
drops rapidly towards higher L. In combination, a flux- 
limited sample of AGN will select many more massive 
BHs in lower mass galaxies than vice versa, biasing the 
sample towards higher Mbh- 

Despite the local comparison sample not being se- 
lected by Mbh - it consists o f inactive galax i es - w e do 
not expect a large b ias: The iMerloni et al.1 (|2009f ) and 
iTrump et al.l (I2009bf) BH masses follow the mass calibra- 
tion bv lOnken et aiT(|2004l ). Their virial M B h for active 
galaxies are calibrated by a forced match to the Mbh-c*- 
relation for inactive galaxies. This might create "wrong" 
Mbh estimates, but it compensates most of the expected 
offset AMbh between our AGN-sample and any (active 
or inactive) sample in the local universe. Together with 
only a small scatter in the local Mbh^M*. bulge it is not 
surprising to find no discernable signature of the LF-bias 
in our data. 

5.2. For Massive Galaxies: Relative Non-Evolution of 
BH Mass Versus Total Stellar Mass 

Hari ng fc Rbd (|2004f ) derive dynamical bulge masses - 
we derive total stellar masses of the host galaxies. Inter- 
preting the coincidence of these two relations at z = 
and z ~ 1.4 depends on how much stellar mass in our 
sample is actually part of a disk component. The stel- 
lar masses of logM^totai = 11.1 ±0.3 are already high in 
the galaxy mass function and our sample galaxies are the 
likely progenitors of giant ellipticals in the local universe. 
If they were already bulge-dominated at the observed 
redshifts, nine billion years earlier, our total masses re- 
vert to bulge masses and our observations mean zero evo- 
lution (though consistent with evolution of factor <2.65 
within errors) in MBH-M^buige, at these A/bh- 

5.3. Evolution of BH Mass Versus Bulge Mass and 
Velocity Dispersion 

While the actual bulge-to-total mass-ratio (B/T) mass 
of our galaxies - in fact of any galaxy at these redshifts 
- is unknown, we have circumstantial evidence that our 
galaxies could contain substantial disk components: (a) 
Visual impression in Figs. QJ+T21 (b) some galaxies have 
Sersic indices near n = 1 (Table [l}, (c) some best fitting 
SSP models have ages down to 1 Gyr (oldest: 5 Gyr), 



(d) the z = 0.36 host galaxies of iTreu et all (120071 ). with 
similar Mbh, have (B/T) mass ~ 1/3, (e) more than 50% 
of log Ma..tnt.a. i = 11.1 ±0.3 galaxi es are strongly starform- 
ing at z ~ 1 ([Noeske et al.ll2007| ) an d (f) of them < 40% 
are on the red sequence (at z = 1.5. lTavlor et~aIll2009T) . 

With a substantial disk component a different inter- 
pretation is possible: 

• We know, independently of their exact evolution- 
ary path, these galaxies have to end up on the local 
Mbh~M*. bulge-relation. If they currently obey the 
same relation but with M*. total, then their masses 
in bulge+disk at z = 1.4 can end up in their bulges 
9 Gyrs later with no addition of mass to either the 
BH or the bulge from outside the galaxy or from 
star-formation. 

• What is required is a process redistributing disk 
stars to the bulge: major and minor galaxy 
merg ing dHopkins et al.l 120091 ) and disk instabili- 
ties ([Parry et alJ 12009) together dominate bulge- 
creation at high masses. 

• In addition to the disk-to-bulge conversion process 
these galaxies can still grow. Merging of similar 
systems will just coadd BHs and stellar components 
separately, moving the system parallel along the lo- 
cal scaling relation. However total mass growth is 
limited since the observed evolution in space den- 
sity of massive ga laxies at z < 1.5 is small (e.g., 
Illbert et~ al. 2009a, and references therein). 

• Even wet mergers are allowed: The gas conversion 
effici ency is not v ery high - in individual merg- 
ers dCrotonl 120061) a nd also overall (at z ~ 0.6, 
iRobaina et alJ I2009D only 10% of SF arise from 
merging. 

• Any subsequent (or ongoing) AGN phase will also 
less than double the BH m ass, since high accre - 
tion states are rare and short ([Hopkins et al.l l2006) . 
Any BH or stellar mass change below a factor of 
two will lie within the 0.3 dex scatter of the local 
relation. 

In the case of a disk component in our galaxies at 
z = 1.4, bulge masses are smaller than the total stellar 
mass in Fig. [3] and the MBH _ M„ s buige-relation will actu- 
ally evolve, when the galaxy-structures are changing over 
time through merging. They will constantly move closer 
to the local relation - consistent with predictio ns from 
simul ations of a merger-driven bulge evolution ([Crotonl 

This has the implication that our result is consistent 
with the non-evolution of the bulge mass relation at 
z < 1.7 for massive early type quasar host galaxies found 
bv lPeng et all ([2006al lbfl - there M* >bu lge = M* )to t a l- At 
the same time we agree with the strong e volution in 
Mrf- M*. b ulge and M-RH~o , *.h iiire claimed by ITreu et al.1 
(2007) and I Woo et alj (|2008T ): If we assume as a limit 
for our galaxies (-B/T) mass < 1/3, our results predict 
an evolution of Mbh/M*, bulge < (1 + z) 1 ' 2 , consistent 
with M BH/M«.bui g c < (1 + z) 15±10 found bv lTreu et all 
(|2007f ) for Seyfert 1 galaxies of the same black hole mass 
at z = 0.36. 
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This path of converging co-evolution of bulge and BH 
is fully indepedent of any interaction or feedback between 
bulge and BH - the AGN can even be switched off since 
z = 1.4. What we are witnessing is the formation of 
the bulge independent of the BH. Only before z — 1.4 is 
a mechanism required to connect total stellar mass and 
BH mass in massive galaxies. We might see an indica- 
tion of this in the large offset of the relat io n for z > 4 
(|Walter et alJ l200i iRiechers et al.l [2008allbl . [2009h . but 
the Lauer-bias and the different mass scale (Mbh > 9.2) 
complicate the picture. Since stellar mass and BH mass 
buildup trace each other very well over cosmic time with 
a factor similar to the local ratio of BH and stellar mass 
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